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A continuous stochastic background of gravitational waves (GWs) for burst sources is produced 
if the mean time interval between the occurrence of bursts is smaller than the average time duration 
of a single burst at the emission, i.e., the so called duty cycle must be greater than one. To evaluate 
the background of GWs produced by an ensemble of sources, during their formation, for example, 
one needs to know the average energy flux emitted during the formation of a single object and the 
formation rate of such objects as well. In many cases the energy flux emitted during an event of 
production of GWs is not known in detail, only characteristic values for the dimensionless amplitude 
and frequencies are known. Here we present a shortcut to calculate stochastic backgrounds of GWs 
produced from cosmological sources. For this approach it is not necessary to know in detail the 
energy flux emitted at each frequency. Knowing the characteristic values for the "lumped" dimen- 
sionless amplitude and frequency we show that it is possible to calculate the stochastic background 
of GWs produced by an ensemble of sources. 
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I. INTRODUCTION 



The detection of gravitational radiation will probably 
mark a new revolution in the history of astronomy. It 
is worth mentioning that the detection of gravitational 
waves (GWs) will directly verify the predictions of the 
general relativity theory concerning the existence or not 
of such waves, as well as other theories of gravity [Q . 

The realm of astrophysics is the place where one finds 
sources of GWs detectable by the GW observatories. 
There is a host of possible astrophysical sources of GWs: 
namely, supernovas, the collapse of a star or star cluster 
to form a black hole, inspiral and coalescence of com- 
pact binaries, the fall of stars and black holes into super- 
massive black holes, rotating neutron stars, ordinary bi- 
nary stars, relics of the big bang, vibrating or colliding of 
monopoles, cosmic strings, etc., among others Nowa- 
days there is a great effort to study, from the theoreti- 
cal point of view, which are the most promising sources 
of GWs to be detected, in particular, their wave forms, 
characteristic frequencies, and the number of sources a 
year that one expects to observe. 

In a few years, instead of building models trying to 
understand how the sources of GWs work, it will be pos- 
sible, starting from the observations (wave forms, ampli- 
tudes, polarizations, etc.) to really understand how the 
GW emission takes place. 

Because of the fact the GWs are produced by a large 
variety of astrophysical sources and cosmological phe- 



nomena it is quite probable that the Universe is pervaded 
by a background of such waves. Binary stars of a variety 
of stars (ordinary, compact, or combinations of them), 
population III stars, phase transitions in the early Uni- 
verse, and cosmic strings are examples of sources that 
could generate such putative background of GWs. 

As the GWs possess a very weak interaction with 
matter passing through it with impunity, relic radiation 
(spectral properties, for example) once detected can pro- 
vide information on the physical conditions from the era 
in which the GWs were produced. In principle it will be 
possible, for example, to get information from the epoch 
when the galaxies and stars started to form and evolve. 

Here we present, in particular, a shortcut to the cal- 
culation of stochastic background of GWs. For this ap- 
proach it is not necessary to know in detail the energy 
flux of the GWs produced in a given burst event. If 
the characteristic values for the dimensionless amplitude 
and frequency are known and the event rate is given it is 
possible to calculate the stochastic background of GWs 
produced by an ensemble of sources of the same kind. 

This paper is organized as follows. In Sec. II we show 
how to calculate the stochastic background of GWs start- 
ing from characteristic values for the dimensionless am- 
plitude and frequency as well as the burst event rate. 
In Sec. Ill we apply the idea presented in Sec. II to 
the calculation of a stochastic background of GWs from 
a cosmological population of black holes, and finally in 
Sec. IV we present the conclusions. 
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II. A SHORTCUT TO THE CALCULATION OF 
STOCHASTIC BACKGROUND OF GWs 

The GWs can be characterized by their dimension- 
less amplitude ft., and frequency v. The spectral en- 
ergy density, the flux of GWs, received on Earth, i^^, 
in erg cm~^s~^Hz~^, is (see, e.g., Refs. 
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where Wobs = 27r!/obs with Vohs the GW frequency ob- 
served on Earth (in Hz), c is the speed of light, G is the 
gravitational constant, and is the strain amplitude 
of the GW (in m'^/'^). For w > 0, Eq.(l) must be mul- 
tiplied by a factor of 2 in order to account for the folding 
of negative frequencies into positive (see, e.g., Rcf. [Q). 
The stochastic background produced by an ensemble of 
sources, of the same kind, would have a spectral density 
of the flux and strain amplitude also related to the above 
equation. The strain amplitude at a given frequency at 
the present time could be, for example, a contribution of 
sources of the same kind but with different masses pro- 
ducing GWs at different redshifts. Thus, the ensemble 
of sources produces a background whose characteristic 
amplitude at the present time is ^/Sh- 

On the other hand, the spectral density of the flux can 
be written as (see, e.g., Ref. @,^) 



(2) 



where fu{voha) is the energy flux per unit of frequency (in 
erg cm~^ Hz~^) produced by a unique source and dR is 
the differential rate of production of GWs by the source. 

The energy flux per unit frequency /ly(fobs) can be 
written as follows (see, e.g., Rcf. 0) 
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With the above equations one finds, for example, the 
dimensionless amplitude of the GWs produced by an en- 
semble of sources of the same kind that generates a signal 
observed at frequency I'obs- Note that in this formulation 
it is not necessary to know in detail the energy flux of 
GWs at each frequency. Knowing the characteristic am- 
plitude for a given source, ft-singio, associated to an event 
burst of GWs, and the rate of production of GWs, it is 
possible to obtain the stochastic background of an en- 
semble of these sources. 

It is worth mentioning that if the collective effect of 
bursts of GWs really form a continuous background the 
quantity called duty cycle must be greater than one. In 
other words: the mean time interval between the occur- 
rence of bursts must be smaller than the typical duration 
of each burst. The duty cycle is defined as follows: 



D{z) = / di?ATGw(l + .z), 



(7) 



where Atgw is the average time duration of single bursts 
at the emission (see, e.g., Ref. Q). 

In the present study we are using the relationship be- 
tween Hbg and hsmgie (Eq. 6) as Ferrari et al. [^^, in 
either case of duty cycle: large or small. In the next sec- 
tion we apply the technique for a case in which the duty 
cycle is small; in another study, to appear elsewhere, we 
apply it to a case where the duty cycle is large (see Ref. 

0)- 

The energy density of GWs is usually written in terms 
of the closure energy density of GWs per logarithmic fre- 
quency interval, which is given by 
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where pc is the critical density {pc — SH"^ /8ttG). The 
above can be written as 
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where /isingic is the dimensionless amplitude produced by 
an event that generates a signal with observed frequency 

l^obs- 

Then, the resulting equation for the spectral density 
of the flux is 



III. APPLICATION: STOCHASTIC 
BACKGROUND OF GWs FROM A 
COSMOLOGICAL POPULATION OF STELLAR 
BLACK HOLES 



From the above equations we obtain for the strain 
1 
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Thus, the dimensionless amplitude reads 
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In this section we apply the formulation presented in 
the precedent section to calculate the background of GWs 
from a cosmological population of stellar black holes. 

From Eq. (6) one sees that it is necessary to know (a) 
ft-singio, here named /ibh, the characteristic amplitude of 
the burst of GWs produced during the black hole forma- 
tion; (b) dR, the differential rate of production of GWs, 
here named (ii?BH , the differential rate of black hole for- 
mation. It is worth noting that we are implicitly assum- 
ing that during the formation of each black hole there is 
a production of a burst of GWs. 
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To proceed it is necessary to know the star formation 
history of the Universe, which we adopted from a study 
performed by Madau et al. which holds for the red- 
shift range < z < 5. It is also necessary to know (a) 
the initial mass function (IMF), which we assume to be 
the Salpeter IMF, and (b) the smallest progenitor mass 
which is expected to lead to black holes (see Refs. |p|,p^). 

A. The rate of stellar black holes formation 

The differential rate of black hole formation can be 
written as 

di?BH = p*{z)--j—(f)(rn)dmdz, (10) 

where /3*(z) is the star formation rate (SFR) density (in 
Mgyr-^Mpc-^), dV is the comoving volume element, 
and (j){m) the IMF (see Refs. 

The SFR density can be derived from observations. In 
particular, our present view of the Universe at redshifts 
z ^ 4 — 5 has been extended by recent data obtained 
with the Hubble Space Telescope (HST) and other large 
telescopes (see e.g. Refs. ||lll-[l^). 

It has been shown that, in general, the measured co- 
moving luminosity density is proportional to the SFR 
density. Thus, the star formation evolution can be de- 
rived from recent UV-optical observations of star forming 
galaxies out to redshifts ^4 — 5 10 . Figure 1 shows the 
SFR density obtained by Madau et al. [|4j. 
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FIG. 1. Evolution of the log of the SFR density 
(Meyr'^Mpc"^) for Qo = I, {A = 0), h = 0.5 and a Salpeter 
IMF. The solid line represents the SFR density evolution 
given by Eq. (11), MDP-1, beyond the dotted line corre- 
sponds to the SFR density given by Eq. (12), MDP-2. 
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FIG. 2. Evolution of the rate of black hole formation occur- 
ring per unit time within the comoving volume out to redshift 
z for Q.0 = 1, (A = 0), h = 0.5 and a Salpeter IMF. The solid 
line represents the rate of black hole formation when we used 
the Eq. (11) to the SFR density, MDP-1, beyond the dotted 
line corresponds to the rate of black hole formation when we 
used the SFR density given by Eq. (12), MDP-2, (see Fig. 
1). 



In particular, there are two different fits to the SFR 
density presented by these authors. The first fit for the 
SFR density (here after referred to as MDP-1) is given 
by 

p*(z) = 0.049 [tl e-*9/0-64 ^ Q 2(1 - e-*^/^-^^)] 

xMq yr"^Mpc"^ (11) 

where tg is the Hubble time in Gyr [tg = 13/(1 -|- z)^/^]. 

The second fit for the SFR density (here after referred 
to as MDP-2) is given by 

/9*(z) = 0.336 e"*«/i-^ -f- 0.0074(1 e-*^/^-^*) 

+0.0197 ^le-'^/O-^"* M© yr^^Mpc"^ (12) 

In the above fits, Eqs. (11) and (12), Madau and col- 
laborators considered an Einstein - de Sitter cosmology 
(r^o = 1) with Hubble constant ffo = 50 kms~^Mpc~^ 
and cosmological constant A = 0. Note that for a dif- 
ferent cosmological scenario it is necessary to rescale the 
SFR density. 

The fit given by Eq. (11) traces the rise, peak, and 
sharp drop of the observed UV emissivity at redshifts 
z ^ 2, while the fit given by Eq. (12) considers that 
half of the present-day stars, the fraction contained in 
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spheroidal systems [Q, were formed at z > 2.5 and were 
enshrouded by dust. This fact produces an increase in 
the SFR density at redshifts z > 2.5 (see Fig. 1) con- 
trary to the sharp drop described in Eq. (11). 

The consistency of p*{z) given by Eq. (12) with the 
Hubble Deep Field (HDF) analysis is obtained assuming 
a dust extinction that increases with redshift. This fact 
is consistent with the evolution of the luminosity density, 
but overpredicts the metal mass density at high redshifts 
as derived from quasar absorbers (see Ref. g). 

Despite this fact, it is interesting also to analyze the 
GW production with p,(z) given by MDP-2 [Eq. (12)] 
because this SFR density produces a large number of su- 
pernovas at z > 2.5, when compared to the SFR history 
described by MDP-1. 

Concerning the IMF here we consider Salpeter's, as 
already mentioned. Thus, 



0(m) = Ato-(i+^\ 



(13) 



where A is a normalization constant and x ~ 1.35 the 
Salpeter exponent. 

The IMF is defined in such a way that (j){m)dm repre- 
sents the number of stars in the mass interval [to, m+dm]. 
The normalization of the IMF is obtained through the re- 
lation 
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where e is the efficiency of generation of GWs and is 
the distance to the source. 

The collapse to a black hole produces a signal with 
frequency (see, e.g., Ref. [Q) 
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where the factor (1 + z)~^ takes into account the redshift 
effect on the emission frequency, that is, a signal emit- 
ted at frequency at redshift z is observed at frequency 

I'obs = 1^0(1 + z)"^- 

From Eqs. (6), (10) and (15) we obtain, for the dimcn- 
sionless amplitude (for a = 0.1), 
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Using this 
obtain A = 



with mi = O.IM0 and m^ = 125Mq. 
normalization of the mass spectrum, we 
0.17 {Mq^-^^. 

In the present work we follow Timmes, Woosley & 
Wheaver ||] (see also Ref. [0), who obtain from stellar 
evolution calculations, the minimal progenitor mass to 
form black holes, namely, ISM© to SOM© depending on 
the iron core masses. Then, we assume that the minimum 
mass able to form a remnant black hole is m^i^ — 25Mq. 
For the remnant mass, M,., we take Mi- = am, where m 
is the mass of the progenitor star and a — 0.1 (see, e.g., 
Refs. [p). 

In Fig. 2 we show the evolution of the rate of black 
hole formation R-Biiiz), i.e., the number of black holes 
formed per unit time within the comoving volume out 
to redshift z, for MDP-1 and MDP-2 for a cosmological 
scenario with SIq = 1-0 and ho — 0.5. Note that MDPl 
and MDP-2 are similar for z < 2.5, and for z > 2.5 they 
are quite different. 



B. The gravitational wave production 

To obtain the stochastic background, besides knowing 
the differential rate of black holes formation presented in 
Sec. II, one needs to know /ibh, the characteristic di- 
mensionless amplitude generated during the black hole 
formation. Following Thorne [|l|, ft-BH reads 



X p,:{z)——(l){m)dmdz 
dz 



(17) 



where in the above equation c^l is the luminosity distance 
to the source. 

The comoving volume element is given by 



Ho J '{l + z)' 
and the comoving distance, r^, is 

_ 2c[l-(l + z)-i/2] 

ITo ■ 



(18) 



(19) 



In the above equation the density parameter is consid- 
ered rio = 1 and Hq is the present value of the Hubble 
parameter. 

The comoving distance is related to the luminosity dis- 
tance by 



dh = ?'z(l + z) 



(20) 



With the above equations we can calculate the di- 
mensionless amplitude produced by an ensemble of black 
holes that generates a signal observed at frequency i-'obs- 

It is worth mentioning that the formulation used here 
is similar to that used by Ferrari et al. 0,^, but instead 
of using an average energy flux taken from Stark & Piran 
[|l6| , who simulated the axisymmetric collapse of a rotat- 
ing polytropic star to a black hole, we use Eq. (15) to 
obtain the energy flux, which takes into account the most 
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relevant quasinormal modes of a rotating black hole and 
represents a kind of average over the rotational param- 
eter. Both formulations present similar results, since in 
the end the most important contributions to the energy 
flux come from the quasinormal modes of the black hole 
formed, which account for most of the gravitational ra- 
diation produced during the collapse process. In a paper 
to appear elsewhere [ |l7| we present a detailed compar- 
ison between our formulation and results with those by 
Ferrari et al. Q. 



C. Numerical results 

Fig. 3 presents the amplitude of GWs as a function 
of the observed frequency obtained from Eq. (17) for 
the two SFR densities present in Sec. II. We obtained 
that the star formation rate given by Eq. (12), MDP- 
2, produces a maximum amplitude /ibg lower than the 
MDP-1 SFR density. This seems to be a contradiction 
since i?BH is higher for MDP-2. Note, however, that for 
z < 2.5, i?BH is higher for MDP-1 and due to this fact 
the maximum amplitude peak is higher for MDP-1. The 
SFR density described by MDP-2 produces a higher i?BH 
for z > 2.5, but the contribution of these events do not 
contribute to enhance the /ibg peak, but instead con- 
tribute to enhance /ibg s-t the lowest frequencies due to 
the redshift effect (see also Fig. 2). 

A comparison of our results with those of Ferrari et al. 
[Ql shows that the formulation used here present similar 
results. It is worth mentioning that in the comparison we 
have adopted e ~ 10~^. Note that in the present study, 
instead of using the average energy flux emitted during 
the axisymmeric collapse of a rotating polytropic star to 
a black hole with different values for the rotational pa- 
rameter, we use Eq. (15). In this equation there is no 
explicit dependence of /ibh on the rotational parameter, 
it represents a kind of characteristic value for the am- 
plitude of GWs during the black hole formation. The 
characteristic frequency given by Eq. (16) has to do with 
the frequency of the lowest m = quasinormal mode of 
a black hole, which is believed to be excited during its 
formation. 

One could argue that it is surprising that our results 
agree so well with those of Ferrari et al., particularly 
those of Fig. 3. The reason for this good agreement is re- 
lated to the fact that the main contribution to the strain 
amplitude in the Ferrari et al. calculations comes from 
the lowest quasinormal mode of the black hole formed, 
which is also the main contribution present in the Eq. 
(15) we use in our study. A close comparison shows, how- 
ever, that the peak of the curve in the Fig. 3 occurs for 
a frequency higher than that of Ferrari et al. Although 
most of the energy comes from the quasinormal mode, 
there is a contribution from the lower frequencies of the 
energy flux (see Rcf. Q|), which moves the peak of the 
strain amplitude of Ferrari et al. to the left as compared 



to ours. 
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FIG. 3. The background amplitude of the GWs as a func- 
tion of Uohs. The solid line (dashed line) represents (sh/e)^^'^ 
for the MDP-1 SFR density (MDP-2 SFR density). 



We show in Fig. 4 ficw as a function of the observed 
frequency. Note that the MDP-1 SFR density presents, 
due to the higher contribution to Rbh for z < 2.5, a 
more relevant contribution to the closure energy density 
of GWs for almost all frequencies, than the MDP-2 SFR 
density does. 

Comparing Fig. 4 with the corresponding figures of 
Ferrari et al., one notes that their curves are broader 
than ours. This occurs due to the fact the closure energy 
density is directly proportional to the energy flux, and 
therefore more sensitive to their frequency dependence. 
The Ferrari et al. energy flux as a function of frequency 
is broader than we use here, this is why their closure en- 
ergy density as function of frequency curves are broader 
than ours. 

Certainly, modifying the exponent x in Eq. (13) to 
X > 1.35, we will obtain a more steeply falling IMF, cor- 
responding to a lower number of massive stars than that 
obtained using Salpeter's IMF. This produces, as a result, 
a lower rate of black holes formation than that obtained 
here, narrowing the curves present in Figs. 3 and 4. How- 
ever, the agreement with the study performed by Ferrari 
et al. will be still good , since that their study present 
the same dependence to the IMF. Thus both results (and 
models) will be modify in the same way. 
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IV. CONCLUSIONS 

Here we present a shortcut to the calculation of 
stochastic background of GWs. For this approach it is 
not necessary to know in detail the energy flux at each 
frequency of the GWs produced in a given burst event, if 
the characteristic values for the "lumped" dimensionless 
amplitude and frequency are known, and the event rate 
is given, it is possible to calculate the stochastic back- 
ground of GWs produced by an ensemble of sources of 
the same kind. 

Since one knows the dominant processes of GWs emis- 
sion one can calculate the stochastic background of an 
ensemble of black holes. We argue that the same holds for 
other processes of GWs production, particularly those in- 
volving cosmological sources, since the number of sources 
could be big enough to produce stochastic backgrounds. 

We apply this formulation to the study of a stochastic 
background of GWs produced during the formation of a 
cosmological population of stellar black holes. We com- 
pare the results obtained here with a study by Ferrari 
et al. [Q, who take into account in their calculations an 
average energy flux for the GWs emitted during the for- 
mation of black holes obtained from simulations by Stark 
& Piran p^ . Our results are in good agreement. 

For most sources of GWs only characteristics values for 
the dimensionless amplitude and frequency are known, if 
these sources are numerous, a stochastic background of 
GWs could be produced. We argue that the formulation 
presented here could be applied to other calculations of 
stochastic backgrounds as well. 
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